Tailoring Phase Space : A Way to Control Hamiltonian Transport 
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We present a method to control transport in Hamiltonian systems. We provide an algorithm - 
based on a perturbation of the original Hamiltonian localized in phase space - to design small control 
terms that are able to create isolated barriers of transport without modifying other parts of phase 
space. We apply this method of localized control to a forced pendulum model and to a system 
describing the motion of charged particles in a model of turbulent electric field. 
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Transport governed by chaotic motion is a main issue 
in nonlinear dynamics. In several contexts, transport in 
Hamiltonian systems leads to undesirable effects. For 
example, chaos in beams of particle accelerators leads to 
a weakening of the beam luminosity. Similar problems 
are encountered in free electron lasers. In magnetically 
confined fusion plasmas, the so called anomalous trans- 
port, which has its microscopic origin in a chaotic trans- 
port of charged particles, represents a challenge to the 
attainment of high performance in fusion devices. One 
way to control transport would be that of reducing or 
suppressing chaos. There exist numerous attempts to 
cope with this problem of controlling chaos 0]. How- 
ever, in many situations, it would be desirable to control 
the transport properties without significantly altering the 
original setup of the system under investigation nor its 
overall chaotic structure. 

In this article, we address the problem of control of 
transport in Hamiltonian systems. We consider the class 
of Hamiltonian systems that can be written in the form 
H = Ho + eV that is an integrable Hamiltonian Hq (with 
action-angle variables) plus a perturbation eV. For these 
Hamiltonians we provide a method to construct a control 
term / of order e 2 with a finite support in phase space, 
such that the controlled Hamiltonian H c = Hq + eV + f 
has isolated invariant tori. For Hamiltonian systems with 
two degrees of freedom, these invariant tori act as bar- 
riers in phase space. For higher dimensional systems, 
the barriers of transport can be formed as a localized 
collection of invariant tori. The idea is to slightly and 
locally modify the perturbation and create isolated bar- 
riers of transport without modifying the dynamics inside 
and outside the neighborhood of the barrier. Further- 
more, we require that, in order to compute the control 



term, only the knowledge of the Hamiltonian inside the 
designed support region of control is needed. 
The main motivations for a localized control are the fol- 
lowing : Very often the control is only conceivable in 
some specific regions of phase space (where the Hamil- 
tonian can be known and/or an implementation is possi- 
ble). Or, there are cases for which it is desirable to stabi- 
lize only a given region of phase space without modifying 
the chaotic regime outside the controlled region. This 
can be used to bound the motion of particles without 
changing the phase space on both sides of the barrier. 

Our algorithm for a localized control contains three 
steps : a global control in the framework of Refs. PlEl- 
a selection of the desired invariant tori to be preserved, 
and a localization of the control term. 
The global control term / of order e 2 we construct is 
such that the controlled Hamiltonian given by H c = 
Ho + eV + f is integrable or close to integrable, i.e. such 
that H c is canonically conjugate to Hq up to some cor- 
rection terms. 

Let us fix a Hamiltonian Hq. We define the linear op- 
erator {H } by {H }H — {H ,H}, where {• ,•} is the 
Poisson bracket. The operator {Hq} is not invertible. We 
consider a pseudo-inverse of {Ho}, denoted by T, satis- 
fying 



{H } 2 r = {H }. 



(1) 



If the operator T exists, it is not unique in general. For 
a given T, we define the resonant operator 1Z as 



71=1- {H }T, 



(2) 



We notice that Eq. QJ becomes {H }TZ = 0. A conse- 
quence is that any element 1ZV is constant under the flow 
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ofH . 

Let us now assume that H is integrable with action- 
angle variables (A, tp) e M z x T' where T' is the l- 
dimensional torus. Moreover, we assume that Hq is lin- 
ear in the action variables, so that Ho = us ■ A, where 
the frequency vector uj is any vector of W . The operator 
{H } acts on V = £ ke » F k (A)e* k as {i/ }^(A, <p) = 
^ k iu; • k Vk(A)e' ik ' ¥ ' . A possible choice of T is then 



rv(A, v )= 



Vk(A) 
iu> • k 



We notice that this choice of T commutes with {Ho}. 

The operator TZ is the projector on the resonant part 
of the perturbation: 



TZV 



Vk(A)e ik -*\ 



(3) 



■k=0 



^From the operators T and TZ, we construct a con- 
trol term for the perturbed Hamiltonian Ho + V , i.e. 
we construct / such that the controlled Hamiltonian 
H c = Ho + V + f is canonically conjugate to Ho + TZV. 
This conjugation is given by the following equation 



i {rv} {H + V + f) = Ho+ TZV, (4) 



where 



/ = E?^V rvr( ^ +1)K (5) 

71=1 

As a consequence, if V is of order e, the largest term in 
the expansion of / is of order e 2 . We notice that Ho and 
TZV are two conserved quantities for the flow defined by 
Ho + TZV (even if H + TZV does depend on the angles in 
general). Therefore, for Hamiltonian systems with two 
degrees of freedom, Ho + TZV as well as the controlled 
Hamiltonian H c are integrable. For higher dimensional 
systems, ifa; is non- resonant, i.e. Vk 6 Z ; \{0}, u> k ^ 0, 
H c is also integrable since it is canonically conjugate to 
Ho + TZV which only depends on the actions. 

Therefore the control term given by Eq. © is able to 
recreate invariant tori, the ones of Hq+TZV. Truncations 
of this control term are also able to recreate invariant 
tori [2|. Of course, the closer to integrability the more 
invariant tori are created. After the computation of the 
control term, the second step is to select a given region 
of phase space where the localized control acts. This 
region has to contain an invariant torus created by the 
previous control and also a small neighborhood of it. 
The invariant torus to be created can be selected by 
its frequency using Frequency Map Analysis 4]. The 
third step is to multiply the control term by a smooth 
window around the selected region where the control 
has to be applied : The locally controlled Hamiltonian is 
H c = Hq + V + x/j where x is a characteristic function 
of a neighborhood P of the selected invariant torus i.e. 




FIG. 1: Poincare surface of section of Hamiltonian © with 
e = 0.034 (enlargement in the inset). 



x( x ) = 1 f° r x £ P, x( x ) 7^ in a small neighborhood of 
P in order to have a smooth function, and x{ x ) = oth- 
erwise. The phase space of the controlled Hamiltonian 
H c is very similar to the phase space of the uncontrolled 
Hamiltonian (since the control term acts only locally) 
with the addition of the selected invariant torus. 
The justification for the localized control follows from 
the KAM theorem : The controlled Hamiltonian can be 
written as H c = H Q + V + f + {\ - 1)/. Since H + V + f 
has the selected invariant torus (by construction), also 
the controlled Hamiltonian H c has this invariant torus 
provided that the hypothesis of the KAM theorem are 
satisfied, namely that the perturbation (x — 1)/ is 
sufficiently small and smooth in the neighborhood of the 
invariant torus. 

The first application of the localized control is on the 
following forced pendulum model Q 



H(p, x, t) = -p 2 + e [cosx + cos(a; — t)] . 



(G) 



for e = 0.034. We notice that for e > 0.02759 there are 
no longer any invariant rotational (KAM) tori [(J. 

For the construction of the control term, we choose 
H = E + cop and V — \fe [p 2 /2 + cos x + cos(x — t)] af- 
ter a translation of the momentum p by tv, and a rescal- 
ing (zoom in phase space) by a quantity \fe. The control 
term of Hamiltonian © becomes 



„ , e 2 /cos x cos(x-t) 

h(x,t)=--( + 1 

2 V ui lu — 1 



(J) 



around a region near p ^ lu. In the numerical imple- 
mentation, we use u> — (3 — -\/5)/2. We notice that the 
perturbation has a norm (defined as the maximum of its 
amplitude) of 6.8 x 10~ 2 whereas the control term has 
a norm of 2.7 x 10~ 3 for e = 0.034. The control term 
is small (about 4%) compared to the perturbation. We 
notice that it is still possible to reduce the amplitude of 
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FIG. 2: Poincare surface of section of Hamiltoniari @ with 
the approximate control term 10 with e = 0.034 (enlargement 
in the inset) using the same initial conditions and same time 
of integration as in Fig. Q 



the control (by a factor larger than 2) and still get an 
invariant torus of the desired frequency. 

The phase space of Hamiltonian © with the approx- 
imate control term {JJ for e = 0.034 shows that a lot of 
invariant tori are created with the addition of the control 
term 0. Using the renormalization-group transforma- 
tion we have checked the existence of the invariant 
torus with frequency u> for the Hamiltonian H + f% with 
e < 0.06965. 

The next step is to localize f 2 given by Eq. around 
a chosen invariant torus created by f% : We assume that 
the controlled Hamiltonian H + f 2 has an invariant torus 
with the frequency uj. We locate this invariant torus us- 
ing Frequency Map Analysis Q . Once the invariant torus 
has been located, we construct an approximation of the 
invariant torus for the Hamiltonian H + /2 of the form 
P = Po(x,t). Then we consider the following localized 
control term : 



A L) = f2(p,x,t)x(p - po(x,t)), 



(8) 



where x is a smooth function with finite support around 
zero. More precisely, we have chosen x(p) = 1 f° r 
H < a , x(p) — for \p\ > P and a polynomial for 
\p\ € [ct,(3] for which x is a C 3 even function. The 
function pg and the parameters a, (3 are determined 
numerically [a — 5 x 10 -3 and /3 = 1.5a). The measure 
of the support of \ is about 1% of the measure of the 
support of the global control. The phase space of the 
controlled Hamiltonian appears to be very similar to 
the one of the uncontrolled Hamiltonian (see Fig. |2J). 
We notice that there is in addition an isolated invariant 
torus which is the one where the control term has been 
localized, i.e. its frequency is equal to (3 — v5)/2. 

The second example comes from plasma physics. It 
aims at modeling chaotic E x B drift motion. In the guid- 



ing center approximation, the equations of motion of a 
charged particle in presence of a strong toroidal magnetic 
field and of a nonstationary electric field are : 



d(x\ c_ ^f-d y V(x,y,t) 



dt \y 



where V is the electric potential, E = — VV, and B = 
Be z . The spatial coordinates x and y where (x,y) 6 K 2 
play the role of the canonically conjugate variables and 
the electric potential V(x, y, t) is the Hamiltonian of the 
problem. A model that reproduces the experimentally 
observed spectrum Q has been proposed in Ref. . We 
consider the following explicit form of the electric poten- 
tial 



A' 



V(x,y,t) = - J2 



1 



x 



r ^i 1 (n 2 +m 2 ) 3 / 2 

n 2 +m 2 <N 2 

sin [2ir(nx + my — t) + ip nrn ] , (9) 

where the phases (p nm are chosen at random in order to 
model a turbulent electric potential. The control term 
has been computed in Ref. [2j from Hq = E [i.e. a local- 
ized control using ui = (1,0)] where E is the conjugate 
action to the angle t mod 2tt : 



8t n f^ {n\ + m\f/ 2 {nl + to 2 ) 3 / 2 

x sin[27r[(n 1 - n 2 )x + (mi - m 2 )y] + (5„ imi - <p n2 m 2 ]- 

(10) 

In Ref. Q , we showed that this global control term is able 
to strongly reduce the chaotic diffusion of test particles 
provided that e < 1, i.e. the trajectories of the controlled 
Hamiltonian V + f 2 diffuse much less than the ones of the 
Hamiltonian V. 

We now localize the control term in a particular region 
of the phase space. We choose this region to be a circu- 
lar annulus of radius ro around a center (a;o, yo) in order 
to determine an area of confinement where the control 
can be applied on the circular border of an experimen- 
tal apparatus. We consider the distance from a point 
(x, y) to the border of the apparatus (circle of radius ro) 
: d(x, y) =| r - y/(x- x ) 2 + (y - y ) 2 \- The localized 
control term is given by 



f2 L) ( x -.y) = h{x,y)x{d{x,y)), 



(11) 



where \{d) = 1 if d < a, x(d) = if d > (3 and x equals 
to a polynomial of degree three for d E [a, 0\ in order 
to ensure the derivability of the controlled Hamiltonian. 
This means that there is a region (the gray circular an- 
nulus in Fig. I3J) where the full control term is applied 
and two border regions (in black) where the control term 
decreases to zero. In Fig. [21 we show some trajectories 
of particles in the cases without control term, with the 
full control term and with the localized control term for 
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FIG. 3: Poincare sections of (a) Hamiltonian @, (6) with 
localized control term ill 1 li and (c) with a full control term 
IIUI I for e = 0.7 using the same initial conditions and same 
time of integration. An enlargement of (b) is shown in (d). 
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FIG. 4: Number of particles Ni„ that remain inside the annu- 
lus during a time t = 30000 by the total number Ntot of parti- 
cles versus the amplitude e of the electric potential © without 
control (open squares), with control term 1101 (open circles), 
and with the localized control term Hit (up-triangles) . 



Hamiltonian with e = 0.7. Whereas in the case with- 
out the control term a lot of particles exit the apparatus 
(the region Q delimited by the circular annulus) and in 
the case with the full control term the trajectories are 
very well confined. The localized control acts only on 
the border without modifying the chaotic motion inside 
and outside the region but catching the particles when 
they arrive close to the border. A measure of this effect is 
provided by the ratio of the number Ni n of particles that 
remain inside the barrier after a given interval of time, 
divided by the total number N to t of particles. In Fig. 0] 
the values of the ratio Ni n /N tot are plotted versus the 
amplitude e of the electric potential J^Jl. It is shown that 
also this localized control is able to significantly reduce 
the escape of particles outside the barrier, although it is 
less efficient than a full control which acts globally on 
the system. We notice that the reduction of transport is 
not achieved by the creation of a KAM torus as it is the 
case for the forced pendulum but by creating a selected 
region of phase space where the system behaves much 
more regularly. 

A measure of the relative size of the control terms 
is given by the electric energies, denoted by £, e 2 and 
associated with the electric potentials V, f% and 
/a , respectively. We define an electric energy £ = 
( | E | 2 ) ■ |fi|/87T where E(x,y,t) = -W and |fi| is 
the area where the control acts. For the potential © 
with N — 25 and for a circular annulus of thickness 
0.7, the relative amplitudes are: ' /£ ~ 0.03 x e 2 and 

e 2 L V e 2 ~ 0.28. For instance, by choosing e = 0.7, 
is 1% of £. This comparison shows that the localized 
control of transport is energetically more efficient than a 
global control of chaos. 

In summary, we have provided a method of control of 
transport in Hamiltonian systems by designing a small 
control term, localized in phase space, that is able to cre- 
ate barriers of transport in situations of chaotic regime. 
Furthermore we remark that in view of practical appli- 
cations, one main feature of our method of localized con- 
trol is that only the partial knowledge of the potential on 
the designed support region is necessary. Our approach 
opens the possibility of controlling the transport prop- 
erties of a physical system without altering the chaotic 
structure of its phase space. 
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